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摘要 : 在 物理 试验 中 准确 测量 控制 棒 价值 对 于 反应 堆 安 全 运行 有 着 非常 重要 的 意义 。 开 展 基于 六 角形 几何 的 多 群 时 
空 动力 学 计算 方法 研究 ， 基 于 此 ， 开 展 适 用 于 六 角形 几何 的 先进 动态 刻 棒 方法 研究 。 该 方法 适用 范围 广 ， 测 量 精度 
高 且 测量 过 程 快 。 利 用 田 湾 VVER 动态 刻 棒 试验 数据 对 该 方法 进行 了 验证 ， 结 果 表 明基 于 该 方法 所 测量 的 控制 棒 价 
值 与 理论 值 偏差 在 10% 以 内 ， 计 算 精度 较 高 ， 且 该 方法 计算 时 间 较 短 ， 计 算 效率 高 。 
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Abstract: Accurate measurement of the control rod worth in physical tests is of great 
significance for the safe operation of the reactor. The research on the calculated method 
of multi~group Space-time neutron kinetics based on hexagonal geometry Was carried out, 
and the advanced dynamic rod worth measurement suitable for hexagonal geometry was studied. 
This method owns high accurate and fast, which has a wide range of applications. The results 
show that the deviation between the control rod worth measured by the method and the 
theoretical value is within 10%, which demonstrated that the calculation accuracy of this 
method is high. The results also show that the calculation time of the method is short, 
and the calculation efficiency is high. 
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控制 棒 价值 测量 是 核电 厂 、 动力 堆 启动 物理 试验 的 一 项 重要 内 容 , 能 很 好 的 验证 堆 芯 性 能 和 设 
计 预 计 值 的 一 致 性 。 动 态 刻 棒 方 法 “是 一 种 快速 的 控制 棒 价 值 ( 以 下 简称 棒 价 值 》 测 量 方法 ， 它 
由 西屋 公司 专家 赵 荣 安博 士 提 出 。 该 方法 最 大 的 优点 是 在 保证 安全 性 的 前 提 下 每 次 缩短 大 修 周期 约 


15 h。 该 方法 目前 已 经 在 中 、 美 、 法 等 国 普及 。 目 前 ， 动 态 刻 棒 研 究 主要 集中 在 基于 方形 组 件 的 堆 


忌 上 ， 对 六 角形 堆 忌 的 应 用 还 较 少 。 

本 文采 用 自主 研发 的 下 XNK 程序 "进行 试验 过 程 的 模拟 ， 计 算得 到 静态 空间 因子 〈SSF) 和 动 
态 空间 因子 (DSF)， 并 且 在 田 湾 VVER 堆 上 完成 了 动态 刻 棒 方 法 的 控制 棒 价 值 测量 试验 。 
1 动态 刻 棒 原 理 

动态 刻 棒 整 个 计算 流程 可 分 为 三 个 步 又 :静态 空间 因子 修正 、 瞬 态 模拟 插 棒 过 程 和 点 堆 逆 动 态 。 
静态 空间 因子 是 为 了 考虑 控制 棒 棒 位 对 组 件 功 率 分 布 和 堆 外 探测 器 的 影响 。 瞬 态 模 拟 揪 棒 过 程 是 为 
了 求 出 揪 棒 过 程 中 探测 器 电流 的 变化 。 点 堆 逆 动 态 是 为 了 求解 控制 棒 下 揪 这 个 过 程 引入 的 反应 性 ， 
称 之 为 动态 反应 性 。 动态 反应 性 不 同 于 静态 反应 性 , 这 是 因为 瞬 态 功率 分 布 并 不 等 同 于 静态 功率 分 
布 ， 即 使 棒 位 一 致 。 在 实际 的 实验 过 程 中 ， 可 以 利用 实验 电流 数据 推出 动态 反应 性 。 动 态 刻 棒 方法 
关键 的 过 程 为 模拟 瞬 态 插 棒 过 程 ， 精 确 地 获得 各 个 节 块 的 功率 变化 。 
1. 1 静态 空间 因子 修正 
首先 建立 探测 器 的 电流 信号 与 堆 蕊 三 维 功率 分 布 的 关系 ， 如 式 (1) 所 示 

R=Cxnx > (WxP xV.) (1) 


其 中 ，R 为 探测 器 信号 ，n 为 堆 芯 功率 ，P 为 节 块 归 一 化 功率 密度 ，W 为 探测 器 响应 因子 ，V 为 节 块 


体积 ，C 为 常数 ， 下 标 k 为 节 块 编号 。 
为 了 考虑 不 同 棒 位 下 功率 分 布 的 差异 ,定义 SSF 为 静态 空间 修正 因子 (Static Spatial Factor) 

> 7 山 
DW x BsnoxV) 


SSF = (2) 


其 中 ARO 表示 全 提 棒 状态 。 

计算 静态 空间 因子 ,得 首先 计算 探测 器 响应 因子 ,通常 采用 MCNP 或 者 确定 论 屏 项 软件 计算 。 然 
后 用 程序 计算 不 同 棒 位 下 的 功率 分 布 ,通过 式 (2) 最 终 得 到 静态 空间 修正 因子 。 注 意 ， 静 态 响应 因 
子 是 棒 组 和 棒 位 的 函数 。 


1.2 动态 空间 修正 因子 
动态 空间 修正 因子 是 逆 动 态 计算 得 到 的 反应 性 与 静态 反应 性 的 比值 ， 即 式 (3)。 


r 
DSF= % (3) 
pa 


为 动态 反应 性 。 动 态 空间 修正 因子 突出 了 稳 态 反应 性 与 动态 反应 性 


其 中 /为 稳 态 反应 性 ，/ 


dyn 


的 差异 。 


2 六 角形 多 群 时 空 动力 学 方程 求解 
2. 1 六 角形 节 块 通 量 展开 

基于 扩散 方程 形式 的 中 子 固定 源 问题 可 以 表示 表示 成 如 下 方程 : 

=Y “D, (r,t ) VD, r,t +E (r,t DPD, (r,t,) 宇 S,(r,t,) 


每 个 节 块 在 径 向 上 都 是 六 角形 , 且 假 定单 个 节 块 内 的 材料 是 均匀 的 。 首先 将 节 块 通过 线性 映射 
E 米 的 参考 节 块 ， 见 图 1。 映 射 由 公式 (5) 表示 : 


(4) 


变换 为 一 个 对 边 距 为 2 厘米 和 高 度 为 2 忆 


式 中 及 有 ,分 别 为 原 节 块 的 对 边 距 和 轴 疝 高 度 。 


图 1 维 参 考 六 角形 节 块 


Fig 1 two-dimensional hexagonal node for reference 


原 固 定 源 问题 (4) 变换 为 


4 0 0 4 0 
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+> (Gt DP, (St)= 5. th $= (05,6,,6:) 
将 每 群 的 中 子 通 量 分 解 为 两 个 部 分 ， 齐 次 项 和 特殊 项 。 齐 次 项 中, 满足 方程 6) 的 齐 次 形式 ， 
I 余弦 和 双 曲 正弦 函数 。 特 殊 项 中 ,用 二 阶 正 交 多 


意义 下 满足 带 源 项 的 方程 (6)。 由 于 每 个 能 群 的 通 量 表达 式 形式 类 似 ， 为 


项 式 展开 ， 并 使 得 在 某 
了 简化 起 见 ， 在 下 面 的 公式 中 如 无 特殊 说 明 ， 将 能 群 下 标 省 略 掉 。 这 样 ， 中 子 通 量 写成 齐 次 项 和 特 


殊 项 的 之 和 ， 见 (7)。 


即 右 端 源 项 为 0， 并 选取 齐 次 方程 的 解 形 式 为 双 


DPD=D,+®, -Deosh( BJ, + Dsinh(Bu + Dwd (7) 


k=1 
as bs di 是 待 求 的 展开 系数 ，VNW 为 齐 次 项 的 数目 ，w 为 正 交 多 项 式 。 
根据 节 块 之 间 的 耦合 条 件 以 及 式 〈6) 确定 展开 系数 。 当 固定 源 求解 结束 后 ， 再 计算 出 先驱 核 
浓度 。 


2.2 向 后 Euler 法 
中 子 时 空 动力 学 方程 可 写 为 如 下 形式 


O 
了 D)—V:D, (r,t)VD, b+2ZebOD。G tb 
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= > (7, DD, (r,t)+ > Yagi G(r,t)+ (1— 3 过 (r,t)D, (r,t) 


(8) 

二 CD= BY 之 nrbg， (r,t)—AC(r,t),i=1,2,...,7 
nl (9) 

将 上 述 两 个 方程 中 的 时 间 导 数 项 采用 向 后 差分 近似 : 

| _ B,D,h 

ot Ab (10) 
OCT _C(r,t,)—C,(r,t,) 

ot ty At a 


离散 得 到 如 下 固定 源 形式 的 方程 组 ， 


1 
—V:D, (r,t )V®D, (r,t )+ 
(r,t VD, (r,t,) 去 


oa =S(r,t,) (12) 
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(13) 


因此 可 以 采用 上 述 介绍 的 节 块 法 求解 。 
2. 3 向 后 微分 公式 


图 2 ”上 半 平 面 的 堆 


法 。 隐 式 线性 多 步 法 


心 布置 图 


Fig 2 The half core layout im the horizontal plane 


一 步 的 计 


向 后 微分 公式 〈BDF) 是 一 类 适合 求解 刚性 问题 的 隐 式 线性 多 步 
算 量 通常 和 向 后 Euler 法 的 计算 量 相当 ， 但 计算 精度 高 于 向 后 Euler 法 。 本 文 考虑 二 阶 BDF 求 解 中 子 
时 空 动力 学 方程 。 对 于 等 步 长 情况 ， 二 阶 BDF 的 计算 公式 由 下 式 给 出 : 
5 一 2 + = At. f,,» (14) 
中 fx 是 了 (biyroy4rt) 的 缩写 ，yix 表示 y(t)。 
3 动态 刻 棒 方法 的 验证 
3.1 程序 验证 
该 基准 题 是 VVER440 堆 型 的 弹 棒 问 题 加 。 堆 芯 材 料 布置 在 水 平 上 关于 x 轴 对 称 ， 见 图 2， 轴 
向 上 的 布置 见 图 3， 其 中 图 中 编号 为 1、2、3 是 燃料 组 件 ， 编 号 为 21、23、25、26 是 控制 棒 组 件 。 
该 算 例 模拟 不 带 反 馈 不 对 称 控制 棒 的 弹出 事故 ， 弹 出 的 控制 棒 价 值 大 约 是 0.78 。 在 初始 时 刻 ， 编 
号 23 和 25 的 控制 棒 并 没 插 入 堆 芯 ， 而 编号 为 26 和 21 的 控制 棒 插 入 堆 芯 ， 棒 位 见 图 3。 初始 时 刻 
编号 26 的 控制 棒 快 速 弹出 堆 芯 ， 且 运动 速度 为 2500cm/s。 在 1 秒 后 ， 堆 芯 紧 急 停 堆 ， 除 了 26 号 
外 的 所 有 控制 棒 〈 包 括 21 号 控制 棒 ) 开始 落 入 堆 芯 底 部 ， 且 以 恒定 的 速度 25cm/s 下 沙 。 整 个 瞬 态 
过 程 持续 6 秒 。 


图 3 扒 芯 在 轴 向 上 的 布置 


Fig 3 the core axial layout 
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图 4 相对 功率 随时 间 的 变化 


Fig 4 relative power varies with time 


3.2 田 湾 VVER 动态 刻 棒 验 证 


根据 田 湾 VVER 三 个 循环 的 实际 


流 数据 ， 利 用 基于 定 步 长 的 BDF 动力 学 模拟 得 到 的 动态 空间 


修正 因子 ， 再 根据 实测 电流 得 到 各 棒 组 价值 ， 实 测 棒 价 值 与 理论 棒 价 值 的 相对 偏差 见 表 1。 从 表 1 


来 看 ， 最 大 棒 组 价值 误差 均 不 超过 10%， 大 部 分 棒 组 价值 符合 很 好 。 
表 1 动态 刻 棒 实 测 棒 价 值 与 理论 棒 价 值 的 相对 偏差 
Table 1 Relative error between the measurement values and the theoretical values of the control rods 
棒 组 第 十 一 循环 第 十 二 循环 第 十 三 循环 
1 9. 41% 6. 2% 1. 7% 
2 —0. 67% 05.0% —1. 4% 
3 0. 50% —2. 4% —9. 2% 
4 —5. 59% 2. 0% 一 7. 4% 
5 —5. 09% —2. 2% —7. 0% 
6 -3. 43% —1. 3% —4. 2% 
7 —3. 48% —3. 0% —4. 1% 
8 —1. 42% 8. 6% 一 2. 4% 
9 0. 95% —1. 3% =2.2% 
10 —3. 2% —4. 7% —4. 0% 
模拟 1 组 棒 下 插 整 个 瞬 态 过 程 ， 在 个 人 PC 机 器 上 ， 平 均 花 费 8 分钟。 这 样 的 计算 效率 允许 整 
个 动态 刻 棒 模 拟 过 程 在 一 个 较 短 的 时 间 内 完成 。 
4 结论 
本 研究 给 出 了 二 阶 向 后 微分 公式 计算 和 函数 展开 粗 网 节 块 法 求解 六 角形 多 群 中 子 时 空 动力 学 
的 理论 ， 并 将 其 与 动态 刻 棒 方 法 结合 起 来 ， 开 发 了 相应 的 动态 刻 棒 程 序 。 通 过 田 湾 VVER 多 个 循环 
的 动态 刻 棒 试 验 数据 表明 基于 所 开发 的 动态 刻 棒 程 序 计算 结果 得 到 的 实测 棒 价 值 与 理论 棒 价 值 的 
误差 均 小 于 10%， 表 明 该 程序 具备 较 高 的 计算 精度 。 另 外 ， 该 方法 的 计算 效率 较 高 ， 计 算 单 组 棒 的 
动态 修正 因子 时 间 约 8 分 钟 。 
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